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ABSTRACT 



Context. HD 50138 presents the B[e] phenomenon, but its nature is not clear yet. This star is known to present spectral variations, 
which have been associated with outbursts and shell phases. 

Aims. We analyze the line profile variability of HD 50138 and its possible origin, which provide possible hints to its evolutionary 
^ stage, so far said to be close to the end of (or slightly beyond) the main sequence. 

. Methods. New high-resolution spectra of HD 50138 obtained with the HERMES spectrograph over several nights (five of them 

'"J ■ consecutively) were analyzed, allowing us to confirm short-term line profile variability. 

' Results. Our new data show short-term variations in the photospheric lines. On the other hand, purely circumstellar lines (such as 

(— I _ [Oi] lines) do not show such rapid variability. The rotational velocity of HD 50138, v^„t = 90.3 ± 4.3 km s"', and the rotation period, 

^ ^ P = 3.64± 1 . 16 d, were derived from the He i /14026 photospheric line. Based on the moment method, we confirm that the origin of this 

I ■ short-term line profile variability is not stellar spots, and it may be caused by pulsations. In addition, we show that macroturbulence 

O ' may affect the profiles of photospheric lines, as is seen for B supergiants. 

^ , Conclusions. The location of HD 50138 at the end of (or slightly beyond) the main sequence, the newly detected presence of line 

C/2 . profile variability resembling pulsating stars, and macroturbulence make this star a fascinating object. 

■ 

I I. Key words, stars: oscillations - stars: mass-loss - stars: winds, outflows - circumstellar matter - stars: individual: HD 50138 

> ■ 

" 1. Introduction tra and discussed several evolutionary scenarios. It seems that 

, HD 50138 might be a star that is close to the end of its main- 

ly . The group of B[e] stars displays, in their optical spectrum, strong sequence evolution. This indicates that this star may represent 
^ . Balmer emission lines and emission from permitted and forbid- ^ jin^ between B[e] and Be stars. Shell phases and outbursts 

■ den transitions of neutral and singly ionized metals. In addition, t^is star have b een reported, based on line profile variabili- 
B[e] stars show circumstellar dust, which is responsible for the jjes (Merrill [79311 Me rrill & Burwell 19331 Hutsemekers [T985l 

" " excess emission in the near and mid-infrared. Especially in the Andrillat & Houziaux IT99T; Pogodin fT997l . And indeed, pho- 
' Galaxy, many B[e] stars are still unclassified with respect to their tometric observations have shown a drop in UBV magnitudes 
evolutionary phase, mainly due to unknown or uncertain stellar ^y about 0.2 mag, which is interpreted as an outburst around 
parameters and especially distances. In addition, they are often 1978/79. In Paper I, it was also suggested that a new shell phase 
embedded in dense circumstellar material. Consequently, their look place prior to 2007 
5, - optical spectra are strongly polluted with circumstellar emis- 

■ sion, while unpolluted photospheric fines, which are needed for ^^^o important to point out that HD 50138 is often dis- 

■ - ■ a proper steUar classification, are rare or absent. '^"^sed in the framework of Herbig Ae/Be stars (see Paper I). 

rru . I IT-, crv no /I r-7,n » J njxir^icox- r .u Bascd On the abscncc of any ncbuloslty around it, as wcU as the 

The star HD 50138 (V743Mon, MWC 158) IS one of these , ^ , . i ■ • • , • 

, . 1 vc J nr T /T 1 irvr,o> 11 1 abscncc of a closc-by star-forming fcgion, it could bc an iso- 

galactic unclassified Biel stars (Lamers et al. 1998), although , ,„ „. i tt , , -^i ■ ■ , 

° rc , , ji j •• 11 lated HAeBiel star However, such a classification is unlikely 

strong ettort has been undertaken to deterrmne Its stellar param- , r.u runu j.i.. uu 

r, . , j ^rvrvr. ! i r r, i ■ bccausc ot thc prcscncc of shell phases and outbursts, which are 

eters. Borges Fernandes et al. ( 12009! hereafter Paper I) reclassi- • i j- , jr,r ^ 

4^ A ,. t)i;TTTT\/ » u A u' u w 1 typical of morc cvolvcd B[e] stars. 

fied It as a B6-7 III- V star based on high-resolution optical spec- ^ 

The presence of a circumstellar disk has been inferred from 

J, _ ' ; '. J ■ , , » T 1 J polarimetric and spectro-polarimetric observations that revealed 

Based on observations made with the Mercator Telescope, operated ^ . . . , f . . . , , • „ 

on the island of La Palma by the Flemish Community at the Spanish ^" i"t"nsic polarization in agreement with a non-sphencally 

Observatorio del Roque de los Muchachos of the Instituto de Astrofi'sica symmetnc distribution of the circumstellar material ( Vaidya 

de Canarias. et al. 1994; Bjorkm an et al. Il998l Oudmaijer & Drew 119991 

** Reduced spectra are only available in electronic form at the Hamngton & Kuhn |2007| l. Recently, Borges Fernandes et al. 

CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via (12011! hereafter Paper II), based on interferometric measure- 

|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/j ments, confirmed the presence of a gaseous and dusty disk 
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around HD 50138. The orientation on the sky plane and a proper 
inclination angle of 56° + 4° were also derived. From both spec- 
troscopic and interferometric observations, the presence of a 
companion is unlikely, although it cannot be definitely excluded 
yet. 

To better understand the nature of its spectral variability, 
we re-observed HD 50138 during several nights (most of which 
were consecutive) in 2011. Thus, we present in Sect. 2 our ob- 
servations and data reduction. In Sect. 3, we present our results 
from these high-quality spectra and in Sect. 4 and 5, we discuss 
the new facts about HD 50138 and present our conclusions, re- 
spectively. 

2. Observations 

High-quality optical spectra were obtained with the High 
Efficiency and Resolution Mercator Echelle Spectrograph 
(HERMES) attached to the 1.2-m Mercator telescope at the 
Roque de los Muchachos Observatory on La Palma (Canary 
Islands, Spain). HERMES is a fiber-fed spectrograph (Raskin et 
al. l201 11 1 and our data were obtained in the high-resolution mode 
(HRF), which provides a spectral resolution of /? ^ 85 000 and 
a field of view of 2.5 arcsec^. The wavelength coverage extends 
from 3770 A to 9000 A. For a proper reduction we used the auto- 
mated data reduction pipeline, HermesDRS, which includes the 
barycentric correction. 

The observations were obtained during eight nights in the 
period 201 1 February 19 to 27. Five of the nights were consecu- 
tive. During each night, for optimum cosmic ray elimination, we 
took two consecutive spectra, which were checked for any line 
variability, and in the case of no change, were then combined 
for a better signal-to-noise ratio {S /N). Details on the resulting 
spectra are summarized in Table 1 . 

Table 1. Journal of available spectra. 



Date 


BJD-2455600 


UT 


T 

' exp 

[s] 


S/N 




[d] 


[h]:[m] 




2011-02-19 


11.5249327 


00:36 


720(2) 


220 


2011-02-20 


12.5581295 


01:24 


900(2) 


140 


2011-02-22 


14.5001431 


00:00 


900(2) 


145 


2011-02-22 


15.4722788 


23:20 


600(2) 


190 


2011-02-23 


16.4782732 


23:29 


750(2) 


240 


2011-02-24 


17.4157888 


21:59 


825(2) 


230 


2011-02-25 


18.4739716 


23:23 


900(2) 


230 


2011-02-27 


20.4782182 


23:29 


1200(2) 220 



Notes. The barycentric Julian date (BJD) and the universal time (UT) 
correspond to the middle of the combined exposures; number in brack- 
ets to the exposure times gives the total number of exposures; signal-to- 
noise (S/N) values refer to the combined spectra. 



3. Results 

3.1. Line profile variability 

The spectrum of HD 50138, as presented in Paper 1, is domi- 
nated by lines from neutral and singly ionized elements. It also 
presents a rich diversity of line profiles from emission to pure 



absorption lines, including several lines with a combination of 
absorption and emission components. 

Line profile variability has previously been reported in the 
literature and associated with shell phases and outbursts (Merrill 
1931; Me rrill & Burwell [19331 Hutsemekers 1985; Andrillat & 
Houziaux ll991t Pogodin[1997]l. Paper 1 reported long-term vari- 
ability in the high-resolution spectra taken eight years apart. The 
variability seen in several lines was interpreted as a new shell 
phase, which took place prior to 2007. 

On the other hand, Pogodin (1997) also reported night-to- 
night variability in the line profiles, from spectra taken in four 
consecutive nights in March 1994, but covering small spectral 
regions. From our HERMES data taken during several consecu- 
tive nights with higher resolution and broader wavelength cov- 
erage, we confirm this short-term variability for a large sample 
of lines of different chemical elements. In Figs.[T]through|9] we 
show some examples of the lines seen in our spectra, noting the 
different strengths of variability, which seem to be related to the 
formation region of these lines. 

The spectrum of HD 50138 displays a few photospheric 
lines. In fact, the only purely photospheric lines seem to be 
He I /14026 and Si ii AAA 128,4131. These lines, as can be seen by 
their residual^ in Figs.[T]and|2l do not display strong variabil- 
ity. These lines are known to be formed at deeper atmospheric 
layers. Other possible photospheric lines like He i A 4009, 4124, 
and 4123 are hardly detected, and HeiA 4388 is blended in our 
spectra. 

On the other hand, lines like the red doublet of 
Si 11/1/16347,6371 (Figs.[3] and |4|, He: lines in the red portion 
of the spectrum, such as He 1/16678 (Fig.|5]) and the core of the 
higher Balmer lines (Figs.|6] and |7]l, present strong variations 
from night to night. These lines are formed in the upper layers of 
the stellar atmosphere or very close to the stellar surface, and are 
extremely sensitive to an enhancement or reduction of the local 
density due to any disturbance in their formation regions. 

The emission lines of circumstellar origin, such as the forbid- 
den lines [Oi] /1/16300, 6364, are quite stable. This can be seen 
in the Fig. [8] for the [Oi] /16364 line and its residuals. Since no 
telluric correction was performed for our data, we do not present 
the [O i] /16300 line, which is polluted by telluric lines. 

However, the same stability is not seen for Ha, as can be 
seen from its residuals (Fig.|9]l. Hq, is usually formed in the outer 
photosphere and lower (ionized) wind or disk region (1-5 7?,), 
while the [O i] lines are formed at much greater distances (> 100 
Rt), where the disk material is already predominantely neutral 
(Kraus et al. l2007ll20T0] l. 

3.2. Rotational velocity 

The stellar rotational velocity projected to the line of sight 
(i> sin /) can be determined from photospheric lines that are nei- 
ther blended nor polluted by emission. As mentioned previously, 
the spectrum of HD 50138 shows only three purely photospheric 
lines: Hei/14026 and Siii/l/l4128,4131. However, the obvious 
large line width results in blending of the Si n lines so that we 
are left with a single line, the He 1/14026 line. It can thus be re- 
garded as the only feasible photospheric line. 

To determine v sin i we applied the well established Fourier 
method (e.g.. Gray [19761 Simon-Di'az & Herrero '2007). The 
Fourier transform of the stellar rotation profile possesses zero 



' All available spectra were added and a mean spectrum was obtained 
and then, the residuals were calculated as the ratio of the individual ones 
with the mean spectrum. 
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Fig. 1. He I /14026 line. The left panel shows its time series and the right panel shows the night-to-night residuals. The time increases 
from bottom (20 11 -02- 19) to top (201 1-02-27), and for a better visual inspection of the spectra, they are slightly shifted along the 
vertical axis with offsets proportional to the interval between the observations, considering the barycentric Julian dates. 






200 400 
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Fig. 2. As in Fig. 1, but for Siii/l/l4128,4131 lines. The radial velocity was derived as a function of the wavelength of rest of the 
A4128 line. 



points, and the first zero is inversely proportional to v sin /. To 
obtain reasonable results, the spectra are required to have both 
a high spectral resolution and a high S IN. While the HERMES 
spectrograph provides an excellent resolution with /? ~ 85 000, 
some of our data suffer from poor S/N. We therefore restricted 
the analysis to four spectra with reasonable S/N values. The 
results are presented in Fig.[TOl The agreement in the position 
of the first zero point in the Fourier transform of these lines 
is excellent. The v sin / values obtained from these zero points 
are Usted in Table|2l and the resulting mean value is usin/ - 
74.7 + 0.8 km s-i. 

Earlier attempts to derive the projected stellar rotational ve- 
locity of HD 50138 have been made by Corporon & Lagrange 
( 1 19991 1 and Fremat et al. (120061 1. who obtained values for v sin / of 
50+30 and 60 + 1 km s ' , respectively. Considering the large er- 
ror bars obtained by these authors, their values are in reasonably 
good agreement with our value. Corporon & Lagrange (il999l l 
measured their value from the He 1/1/14472, 6678 lines, which 
are both highly variable in our spectra (for He 1/16678 line see 
Fig.|5]i. Also Fremat et al. ( 120061 ) used highly variable lines like 



Table 2. Rotational velocities (ysin/) derived from the 
HeU4026 line. 



Date 


V sin ( 


S/N 




[kms-'] 




2011-02-23 


75.39 


240 


2011-02-24 


74.75 


230 


2011-02-25 


73.36 


230 


2011-02-27 


75.30 


220 



the He i and Mg ii lines at AAAAll, 448 1 , and their S /N was con- 
siderably lower than ours with about 150. Our value, which is 
obtained from the He i /14026 line, the stablest line in the whole 
spectrum, can thus be regarded as reliable. 

Considering in addition an inclination of / = 56° ± 4° of the 
star as obtained for its dusty disk from an analysis of interfer- 
ometric observations in Paper II, the stellar rotational velocity 
of HD 50138 becomes v„i = 90.3 ± 4.3 km s^'. Applying the 
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50 100 
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bO 100 150 

Radial Velocity {km s"') 



Fig. 8. As in Fig. 1, but for [O i] ^6364 line. 




-500 500 

Radial Velocity (km s"') 




500 
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Fig. 9. As in Fig. 1, but for Ha /16563 line. 



stellar parameters obtained in Paper I, the rotational period of 
HD50138isP = 3.64+1.16d. 

Next we computed synthetic, rotationally broadened photo- 
spheric line profiles for stellar parameters of Teff = 13 000K 
and log^ - 4.0, in agreement with Paper I, which we fitted to 
the observed He i and Si ii lines, seen in the combined spectrum 
of 201 1-02-27. These are shown as the dashed lines in Fig.fTTI 

Comparing the synthetic and observed profiles, it is obvious 
that stellar rotation alone cannot account for the observed line 
width. This is particularly evident for the Si ii lines (top panel of 
Fig.fTTTi. To fit the observed profiles it was necessary to broaden 
the lines with an additional, Gaussian shaped profile. The needed 
amount of the extra velocity broadening was about 30-40 km s"' . 
Line profile fits using these two values are shown in Fig.fTTI as 
dotted and solid lines, respectively. 

The presence of such an excess broadening was also reported 
for another group of stars, the B supergiants (e.g., Ryans et 
al.|2002; Simon-Diaz & Herrero l2007l Markova & Puis 2003. 
These stars are furthermore known to display strong line profile 
variabiHty (e.g., Kaufer et al. jgg ?, 200 6, Lefever et al. 20^ 
Markova et al. |2008 ; Clark et al. 1201 Oi l, and recent studies by 
Simon-Diaz et al. (.2010i have revealed that a correlation exists 




201 1-02-27 

2011-02-25 

201 1-02-24 

2011-02-23 



Fig. 10. Position of the first zero point in the Fourier transform 
of the line profiles of the He i A4026 line. 



between line profile variability in B supergiants and their excess 
broadening, which is usually referred to as 'macroturbulence'. 
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- [A] 




* [A] 



Fig. 11. Fit of the observed Si ii (top) and He i (bottom) lines with 
synthetic photospheric lines, using stellar parameters in agree- 
ment with Paper I, broadened with stellar rotation and macrotur- 
bulence. 



3.3. Moment analysis 

While the origin of the macroturbulence velocity seen in the 
line widths of B supergiants was long unclear, recent studies 
by Aerts et al. (I2009I I suggest that the broadening might result 
from many (hundreds) individual contributions provided by stel- 
lar pulsations. If pulsations could be a plausible explanation for 
the excess broadening seen in the He i line of the non-supergiant 
B[e] star HD 50138, then this line should also show profile vari- 
ability to at least some extent. Since the residuals of this line do 
not reveal any obvious variability, a more thorough analysis of 
its profile is needed to unveil any small-scale (periodical) vari- 
abiUty. 

To investigate the putatively stable He i line in more detail, 
we make use of the moment method, which is an excellent tool 
for studying even small-scale line profile variabilities and distin- 
guishing between different physical mechanisms as the origin of 
the variabilities, such as (non-)radial pulsations and stellar spots. 
But for a reliable distinction between the two scenarios (pul- 
sations versus chemical surface inhomogeneities in the form of 
spots), the variability in different chemical elements needs to be 
investigated. Besides the He i A4Q26 line, we therefore included 
the Si 11/1/14128, 41 31 lines in the analysis, but keeping in mind 
that the results of these lines probably have a higher uncertainty 
because of their blending. We followed the description of Aerts 
et al. ( I1992I I and North & Paltani ( 1994) by computing the equiv- 
alent width, together with the first three moments of these lines. 



Due to their poor quality, the spectra taken on 2011-02-20 and 
201 1-02-21 were excluded. 

Figure[T2] shows the resulting time variation of the equiva- 
lent widths and the first three moments, which represent the ra- 
dial velocity (< y' >), a measure of the width (< >), and 
the skewness (< >) of the line. Quite clearly, all three lines 
show variability in their first three moments, while their equiva- 
lent widths (zeroeth moment) remain relatively stable. The am- 
plitude in radial velocity variation is on the order of lOkms 
It is slightly higher in the Sin lines. This might explain why 
the residuals in those lines show some signal (see Fig.|2]l, while 
those of the He 1/14026 line are almost featureless. 



4. Discussion 

The very similar variation in both He i and Si ii lines makes it 
unlikely that their variabilities could originate in a star with a 
heterogeneous surface abundance pattern in the form of stellar 
spots. Such stars usually show a quite different behavior in the 
variability seen in their helium lines with respect to those in their 
metal lines (see, e.g.. Briquet et al. 120011120041 Lehmann et al. 
I2OO6I 1. 

In addition, the line profiles in spotted stars vary with their 
rotational period. Owing to an insufficient time coverage of our 
observations (only one data point per night), a search for period- 
icity in the moments failed. However, the rapid change in radial 
velocity from -5 to about +10kms ' in the Hei line within less 
than 26 hours (from BJD 17.416 to 18.473), while it showed a 
steady decrease from ~ 6 to ~ -5kms ' during the three pre- 
ceding nights, might indicate (using a tentative sine curve fit) a 
possible (maximum) period of ~ 1.7 days, but the real period 
could also be (much) shorter. This needs to be studied on time 
series of spectroscopic observations taken with a much better 
time resolution. 

Consequently, the variability in the lines of HD 50138 is also 
not connected with the rotation period of ~ 3.6 days. Instead, it 
resembles those seen in pulsating stars. In addition, the strong 
short-term variability of the line profiles is mainly contained in 
the lines formed in the upper layers of the stellar atmosphere 
or very close to the stellar surface (like the core of the higher 
Balmer lines), while the lines formed in the circumstellar en- 
vironment (like the [O i] lines) are not affected by these night- 
to-night variations. It is interesting to point out that two pulsat- 
ing late-type Be stars exist and occupy a similar region in the 
Hertzsprung-Russell diagram as our object (see Fig.[T3]). 

In addition, the fact that macroturbulence is present in the 
photospheric lines of HD 50138, a known property of B super- 
giants, and its classification as a star at the end of (or slightly 
beyond) the main sequence, makes it a fascinating object. 



5. Conclusions 

Based on new high-resolution spectra of HD 50138 obtained 
during eight nights (five of them consecutive), we confirm the 
presence of short-term variability in a large sample of lines 
of different elements. These spectral variations have different 
strengths depending on the line formation region. The strongest 
short-term variability is identified in the profiles of lines arising 
in the upper layers of the stellar atmosphere or very close to the 
stellar surface. 

Analyzing the few present photospheric lines, we derived the 
rotational velocity and rotational period of HD 50138, Oiot = 
90.3 + 4.3 km s"' and P = 3.64 + 1.16d, respectively. We 
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He I (4026 A) 



Si II (4128 . 




0.5 

0.4 

0.3 

0.2 

0.1 


10 



1500 - 



15 

BJD [2455600+] 



Si II (4131 A) 



15 

BJD [2455600+ 



Fig. 12. Variation in the equivalent widths and first three moments calculated for the He i A4026 and Si ii AAA 128,4131 with time. 
The moments < u' >, < >, and < >, have units (kms"'), (kms"')^, and (kms"')-', respectively. 



conclude that the short-term line profile variability is not con- 
nected with this rotation period. In addition, we used the mo- 
ment method to study the origin of the line profile variability 
and discard any possibility that there are stellar spots, which sug- 
gest pulsations as the origin of the short-term variations identi- 
fied. However, from our data it was not possible to determine the 
modes of these pulsations. Nevertheless, the scenario of pulsa- 
tions, though putative at that moment, might also favor the ap- 
pearance of shell phases as reported in Paper I. As in pulsating 
Be stars, these mass ejections could be associated to possible 
beat periods. 

These new results, accompanied by the probable macrotur- 
bulence in the photospheric lines, reveal the fascinating and 
complex properties of HD 50138, an object clearly in need of 
more detailed study. An observational campaign to obtain more 
data of higher quality, in a longer temporal series, is definitely 
necessary. From these data, it will be possible to precisely de- 
termine the real origin of the line profile variations and conse- 
quently the nature of this star 
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Fig. 3. As in Fig. 1 , but for Si n ^6347 line. 
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Fig. 4. As in Fig. 1 , but for Si ii /1637 1 line. 
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Fig. 5. As in Fig. 1 , but for He i A6678 Une. 
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Fig. 6. As in Fig. 1, but for H 10 ^798 line. 
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Fig. 7. As in Fig. 1, but for H 9 A3835 Une. 



